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Abstract. We present first results from contemporaneous observations of Cygnus X-1 with INTEGRAL and RXTE, made 
during INTEGRAL'S performance verification phase in 2002 November and December. Consistent with earlier results, the 
3-250 keV data are well described by Comptonization spectra from a Compton corona with a temperature of kT ~ 50-90 keV 
and an optical depth of t ~ 1.0-1.3 plus reflection from a cold or mildly ionized slab with a covering factor of Q./2n ~ 0.2- 
0.3. A soft excess below lOkeV, interpreted as emission from the accretion disk, is seen to decrease during the 1.5 months 
spanned by our observations. Our results indicate a remarkable consistency among the independently calibrated detectors, with 
the remaining issues being mainly related to the flux calibration of INTEGRAL. 
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• -h 1. Introduction 

$H The broad band X-ray continuum of Galactic black holes (BH) 
. i in the hard state is characterized by a power-law with photon 
index F ~ 1.7, and an exponential cutoff at ~150keV, which 
is gen erally attributed to Comptonization JSunvaev & Trum pet 
Il979l) . A soft excess below ~1 keV, a hardening above lOkeV 
attributed to Compton reflection, and a fluorescent 6.4 keV Fe 
Ka line are all taken as evidence for the presence of a (possi- 
bly mildly ionized) accretion disk. Long-term monitoring pro- 
grams and the existence of state changes from the hard state 
to the soft state (see, e.g., Remillard 2001) show that the rela- 
tive contribution of the Comptonizing medium and the accre- 
tion disk to the emitted X-ray luminosity are variable on long 
timescales dPottschmidt et alJ2003UZdziarski et all2002l) . 



The knowledge of the spectral shape above lOkeV is cru- 
cial for the determination of the temperature and optical depth 
of the Comptonizing medium, as well as for the determination 
of the reflection fraction. With the launch of the International 
Gamma-Ray Astrophysics Laboratory (INTEGRAL), a new 
instrument providing data from 5 keV to several MeV has be- 
come available. To test the INTEGRAL instruments and fur- 
ther our knowledge of BH spectra, we organized simultaneous 
measurements of the canonical BH Cygnus X-1 with the Rossi 
X-ray Timing Explorer (RXTE) during INTEGRAL 's perfor- 
mance verification (PV) phase. Here, we present first results of 
these observations. In Sect.|5]we describe the observations and 
the data extraction, followed in Sect. [5] by the results of our 
modeling. Sect.|4]summarizes our main points. 



* Based on observations with INTEGRAL, an ESA project with 
instruments and science data centre funded by ESA member states 
(especially the PI countries: Denmark, France, Germany, Italy, 
Switzerland, Spain), Czech Republic and Poland, and with the par- 
ticipation of Russia and the USA. 

Correspondence to: K. Pottschmidt (katja.pottschmidt@obs.unige.ch) 



2. Observations and data analysis 

Cyg X-1 was observed during the PV phase of INTEGRAL 
in a variety of different observing modes. RXTE observing 
times were chosen to coincide with INTEGRAL observations 
that pointed at the source (satellite revolutions 11, 14, 16, 
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Table 1. Log of the simultaneous INTEGRAL /RXTE obser- 
vations, showing the INTEGRAL revolution (Rev.), the date of 
the observation, as well as the exposure time of the instruments 
used in the data analysis in ksec (for RXTE, we give the PCA 
exposure and the dead-time corrected HEXTE live-time). 

JEM 



Rev. 


Date 


RXTE 


XI 


X2 


ISGRI 


SPI 


11 


2002-11-16 


12.1/8.1 






11.7 


86.0 


14 


2002-11-26 


6.7/4.1 




15.0 


0.6 




16 


2002-12-02 


4.1/2.6 




10.1 


5.8 




18 


2002-12-06 


7.5/5.2 


10.6 




12.2 




25 


2002-12-29 


29.7/34.9 






8.0 


86.4 



and 18) as well as with observations during revolution 25, 
where Cyg X-l was on average ~10° off-axis. During revo- 
lutions 14, 16, and 18 INTEGRAL performed staring observa- 
tions dedicated to the calibra tion of the Joint European X-ray 
Monitor (JEM-X; iLund et alJl2003l) and the Imager on-board 
INTEGRAL QBIS: IUbertini et alJ2003l) . For both instruments, 
data simultaneous to RXTE were extracted on 'science win- 
dow timescales' (i.e., stable INTEGRAL pointings of ~2000s 
length). For these revolutions, no source spectrum could be 
extracted from the Sp ectrometer on-board INTEGRAL (SPI; 
IVedrenne etai1l2003l) . as INTEGRAL was not in one of its 
dithering modes and the SPI imaging solution is underdeter- 
mined. Information from SPI is only available for the hexag- 
onal dithering mode in revolutio n 1 1 and for the 5 x 5 dither 

I ~ II 1 

pattern of revolution 25 (see Courvoisier et al. 2003 for a de- 
scription of INTEGRAL'S dithering patterns). Due to the small 
effective area of SPI, we use data from the whole INTEGRAL 
revolution contemporaneous to the RXTE observations. For the 
IBIS analysis of these two revolutions, only the closest con- 
temporaneous science windows were considered. These were 
either on-axis (rev. 11) or had an offset of < 5° (rev. 25). The 
JEM-X data extraction failed for these two revolutions. Tabled 
provides a log of al l th e observations analy zed herein. See 
Bazza no et al and lLaurent et al.l J2003h for discussions 

of (non-simultaneous with RXTE) IBIS data of Cyg X- 1 . 

The INTEGRAL data were extracted using the 
INTEGRAL off-line scientific analysis (OSA) software, 
version 1.1, including those fixes to the released software 
which were available up to 2003 June 15. For all instruments, 
the newest available response matrices were used. Data from 
JEM-X were considered from 10^-OkeV only, due to on-board 
event selection effects in the early phase of the mission (before 
rev. 45). For IBIS we used data from 40-250 keV, the energy 
band in which the energy response of IBIS was considered 
calibrated at the time of writing, and only considered the 
INTEGRAL Soft Gamma-Ray Imager (ISGRI) part of the 
detector. We applied the 10% systematic uncertainty recom- 
mended by the IBIS team (Ubertini, priv. comm.) in quadrature 
to the data. For SPI, data f rom 50-200 keV were ex tracted 
with SPIROS, version 4.3.4 JSkinner & Connellll2003l IStrond 
120031) . See Sect. [3] for a discussion of this rather low upper 
energy threshold of SPI. For all instruments, background 
subtraction has been taken into account by the OSA. Current 
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Fig. 1. a 15 GHz radio lightcurve from the Ryle telescope and b 
2-12 keV RXTE-A11 Sky Monitor flux of Cyg X-l for the time 
of the pointed INTEGRAL and RXTE observations (dashed 
lines). Note the end of the 2002 soft state of Cyg X-l in 2002 
mid-October. 



experience shows, e.g., that with ISGRI a few ks are sufficient 
to extract a significant spectrum up to ~ 1 00 ke V for a 60 mCrab 
source lPaizisetalJ E()03). thus background influences should 
not be an issue for Cyg X-l up to a few lOOkeV. 

We compare the INTEGRAL data with measurements from 
both RXTE pointed instruments , the Proportional Counter 
Array (PCA: Ijahoda et alJ 1 19971) and the High Energ y X- 
ray Timing Experiment ("HEXTE: iRothschild et al!l99 8). The 
RXTE data were extracted with the HEASOFT s oftware, ver- 
sion 5.2, using our standard procedures (e.g., I Wilms et alJ 
Il999h . The high effective area of the PCA (-4000 cm 2 ) ne- 
cessitates the use of energy dependent systematic uncer- 
tainties of 1% for PCA channels 0-15 (<7.5keV), 0.5% 
for channels 16-39 (7.5ke V-18keV), and 2% up to 25keV 
jRrevkenbohm et afll2003l) . These uncertainties are added in 
quadrature to the data. For the HEXTE, the 20keV-200keV 
energy range was considered. Finally, the Ryle telescope in 
Cambridge, England, performed simultaneous radio observa- 
tions at 15 GHz in addition to its daily monitoring of the source 

(Fig.m. 

Spectral fi tting was performed with XSPEC 11.2.0bc 
( Arnaud 1996), which includes important fixes to the reflection 
model used in the spectral analysis (Sect. [3ji. The uncertainty 
of the instrument flux calibration was taken into account by 
introducing a multiplicative constant into the spectral models, 
and normalizing all fluxes to the PCA. This constant was 0.80- 
0.84 for the HEXTE (consistent with earlier results), 0.56 for 
JEM-X 1, 0.43 and 0.46 for JEM-X2, 1.03 and 1.05 for SPI, and 
it varied between 0.89 and 1.38 for ISGRI (see Tables |2]and|3] 
below). The latter large variation is mainly due to known issues 
with the current ISGRI deadtime correction. 
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Fig. 2. Spectra, best fit compTT model, and residuals for the observations of 2002 December 02 and 2002 December 29 (PCA, 
3-25 keV, light blue; JEM-X, lO^lOkeV, green; HEXTE, 20-200 keV, red; IBIS, 40-250 keV, dark blue; SPI, 50-200 keV, olive 
green). For SPI, the spectrum as provided by the SPI spectral extractor, SPIROS, is shown, for the other instruments, we display 
the count rate normalized by the detector effective area. The PCA residuals for revolution 25 are consistent with known systematic 
effects. 
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Fig. 3. Unfolded X-ray spectra for the observations of Fig. |5J and the components of the compTT model. For all instruments, the 
measured fluxes have been renormalized to that measured by the PCA, the colors are those of Fig.|2] 



3. Modeling the 3-250 keV spectrum 



We modeled the spectra using the canonical model for the hard 
state, a Comptonization continuum which is partly reflected off 
a slab-like accretion disk. In these models, it is assumed that 
the Comptonizing medium sandwiches the accretion disk, al- 
though its covering factor can be less than unity. This geometry 
is not the only possible one for Cyg X-l. Sphere+disk mod- 
els, e.g., where the Comptonizing medium is an inner spherical 
accretion flow irradiated by an outer accretion disk (possibly 
slightly penetr ating the sphere) have also been shown to de - 
scribe the data JDove et alJll 99711 1998llZdziarski etaill 19991) . 
Furthermore, emission from th e base of a jet might also con- 
tribute to the observed X-rays jMarkoff etal]l2002l) . We will 
consider such more complex models in a future paper. 



Due to the limitations intrinsic to the available 
Comptonization models, it is generally useful to apply 
several different Comptonization models to the data. Such an 
approach allows one to test whether any trends seen in the 
evolution of spectral fit parame t ers are stable against these lim- 
itations. See, e.g.. lNowak et alJ J200 2) for a discussion of these 
issues in the context of RXTE observations of GX 339-4. We 
therefore describe our joint INTEGRAL-RXTE spectra using 
two different k inds of spectral models . The first one is based on 
the theory of lHua & Titarchukllll995i XSPEC model compTT) 
as a Comptonization continuum, which assumes a Wien 
spectrum for the seed photons. This continuum is then partly 
reflected off a cold accretion disk jMagdziarz & Zdziarski 
l!995t) . To account for the effects of a patchy Comptonizing 
plasm a, an additional soft disk black body ( Mi tsuda et alJ 
1984) is added. The temperature at the inner edge of the 
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Table 2. Best fit parameters for the compTT model (full model in XSPEC notation: 
constxphabs [diskbb+gauss+corapTT+re£lect(comptt)]). The parameters shown are the inner accretion disk tem- 
perature kT m and its normalization, A^k = ((^in/km)/(D/10kpc)) 2 cos i, the electron temperature of the Comptonizing plasma, 
kT e and its optical depth r, the covering factor of the cold reflecting medium, Q/27T, and the energy and width of the Gaussian 
Fe Ka line, E Ka and cr Ka , and the flux normalization constants of the individual instruments with respect to the PCA are given 
by chexte, ciSGRb Oem-xi, cjem-x2, and cspi- For revolutions 11 and 25 the results are shown with and without taking SPI into 
account in order to illustrate its systematic influences. 
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Table 3. Best fit parameters for the eqpair fits (full model in XSPEC notation: constxphabs [diskbb+gauss+eqpair] ). 
Shown are the disk parameters kT m and Aa^, inner accretion disk temperature, kT m , the ratio of the compactness of the 
Comptonizing plasma to the accretion disk compactness, lh/k, from which the Comptonization temperature kT e can be de- 
rived (as kT e is not a fit parameter, we do not give its uncertainty), the optical depth of the Comptonizing plasma, r, the covering 
factor of the reflecting medium, Q./2tt, and its ionization parameter, the energy E^a and width <x of the Fe Ka line (where Em a 
was limited to >6 keV in the fits), and the flux normalization constants of the instruments with respect to the PCA. As in Tab. [2] 
for revolutions 1 1 and 25 the results are shown with and without taking SPI into account. 
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disk, kTi n , is set equal to the seed photon temperature. As the 
reflection model does not include fluorescent line emission, 
the Fe Ka line at 6.4 keV is described by a Gaussian. We 
also use a fixed absorbing colu mn of Nn = 6 x 10 21 cm -2 
( Bal ucihska-Church et aD Il995l) . as the PCA 3keV lower 
threshold does not allow such low columns to be constrained. 



We also described our data using the thermal 
Comptonization model of IConnil ill 9991 XSPEC modelname 
eqpair) 1 , which includes reflection off a possibly ionized 



1 The eqpair-model is available as a local XSPEC model at 
http : //www . astro . vale . edu/coppi/. 
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accre tion disk ( Mag dziarz & Zdziarskil 1 19951 iDone et all 
Il992l) . In this case, the spectral shape of the seed photons is a 
disk black body, part of which is not Comptonized to take into 
account the possibility of a patchy Compton corona. This has 
again been realized in form of an additional disk black body 
component of the same temperature as the seed photons. The 
Fe Kor line is again described by a Gaussian and absorption is 
treated identical as in the compTT fits. 

Both models resulted in equally good descriptions of the 
data. We show examples of our best fits in Figs.|2]and[3l and list 
the best fit parameters for all observations in Tables[5]and[3] For 
revolution 25,x^ ed = 2.1, for the other observations^^ < 1.5. 
Due to the long PCA exposure of rev. 25, x 2 is dominated by 
the PCA (the PCA residuals shown for revolution 25 in Fig. [3 
are consistent with calibration uncertainties). We therefore con- 
clude that all of our Comptonization fits resulted in a satisfac- 
tory description of the 3-200 keV data from Cyg X-l. 

For the compTT model, the Comptonization parameters are 
in principle consistent with earlier observations JDove et alJ 
1998; Gierlihski et alJ Il997h . For aU five observations, the 
Comptonization temperature lies between 50 and 60keV and 
the Thomson optical depth of the corona between 1 and 1.2 
(cf. the two last columns of Table [2] for systematic influences 
from the SPI spectrum). These parameters vary only slightly 
with time, as is expected for the hard state of Cyg X-l. The 
covering factor of the reflecting slab is Q./2n ~ 0.18. This 
rather small value could imply that the Compton corona only 
covers part of the accretion disk, or that the accretion ge- 
ometry is more complicated than the slab geometry assumed 
here ( sphere+disk models, e.g., imply £l/2n < 0.3: lDove et alJ 
119971) . The eqpair models agree qualitatively with these re- 
sults and also require moderate t and Q./2t: < 0.5. They also 
allow for a very modest amount of ionization, with ionization 
parameters £ < 100. For both, compTT and eqpair, iron lines 
were consistent with equivalent widths of ~50 to lOOeV, and 
are possibly broad. 

Contrary to this constant behavior above lOkeV, the spec- 
trum is m uch more variable on the long term at lower ener- 
gies (e.g., IZdziarski et alll2002l Fig. 12). Here we see a sys- 
tematic evolution of the soft X-ray spectrum. For the compTT 
model fits, the flux of the non-Comptonized disk component 
decreases by a factor of about two between the first and the last 
of our observations, whereas for the eqpair fits, the flux of 
the additional disk component nearly completely vanishes. It is 
tempting to interpret this decrease in the flux of the soft model 
component as a true physical indication of the final stages of 
the 2002 transition out of the soft state, given the proximity of 
our first observation to the end of this state (Fig.Q. 

The temperature of the non-Comptonized disk component 
is rather high for the compTT fits (kT m ~ 700 eV, implying an 
inner disk radius comparable to that of the last marginally sta- 
ble orbit of a disk around a 10 M Schwarzschild black hole). 
There have been reports of a soft excess in the hard state of 
Cyg X-l in addition to th e black body provided by the disk 
(e.g. JPi Salvo et aTll200lh and the hard Comptonized compo- 
nent. In principle, such an excess might explain our compTT 
parameters, however, our data do not allow to evaluate more 
complex models for the soft energy regime. Also, the disk tem- 



perature found in the eqpair fits is significantly lower. Here 
it is close to the kT m = 200 eV found, e.g., with ROSAT 
( Bal ucihska-Church et alJll995l) . with an implied disk radius 
near 50 GM/c 2 for a 10M o BH. The Fe Ka line energy in the 
eqpair fits is systematically too low, however, possibly indi- 
cating the soft spectral complexity noted above. 

This discrepancy of the temperatures and implied disk radii 
might be caused by slight differences in the treatment of the 
transition of the seed photon spectrum to the Comptonization 
spectrum, especially given that compTT describes the seed pho- 
tons with a simple Wien spectrum. Note also that constrain- 
ing the temperature of the disk black body component to the 
seed photon temperature - a rather established simplification 
- removes some degeneracy from the low energy parameters. 
Finally, we also cannot exclude the possibility that part of the 
fitted black body is due to the calibrations of the PCA below 
~ 10 keV, especially given that the temperatures of the soft com- 
ponents for both cases, compTT and eqpair fits, imply the peak 
emission is falling below the ~3keV lower cutoff of the PCA. 
We conclude that our simplifications of the low energy spectral 
model are justified, especially since their effect on the param- 
eters derived for the spectrum above lOkeV can be expected 
to be negligible. Note that the spectral shapes measured above 
lOkeV band with the PCA, HEXTE, JEM-X, and ISGRI are 
indeed remarkably consistent with each other, suggesting that 
their respective calibrations are in agreement. 

Measurements with the Compton Gamma-Ray Observatory 
(CGRO) and RXTE show the presence of a hard tail 
from non-thermal Comptonization at >300keV in the hard 
state of Cyg X-l and above several keV in the soft state 
dGierlinski et a lJ Il 999t iMcConnell et al Jl200ol l2002fc . In prin- 
ciple, ISGRI and SPI detect Cyg X-l to at least 500 keV and 
should be able to detect significantly this hard tail. Currently, 
however, the intercalibration of the instruments above 250 keV 
is not yet good enough to allow an extrapolation to this energy 
band. While ISGRI and HEXTE are in good agreement, show- 
ing a spectral turnover and having a peak of v f v at ~100keV, 
the SPI spectrum is consistent with a continuation of a pure 
power-law. 

4. Conclusions 

We have presented first results from fitting simultaneous 
INTEGRAL-RXTE data of Cyg X-l with two different 
Comptonization models. Our results for the spectral shape are 
in basic agreement with earlier results, with HEXTE, ISGRI, 
JEM-X, and to a lesser extent SPI, agreeing with each other. 
We take this result to mean that the energy redistribution of the 
INTEGRAL instruments is well understood up to ~250keV. 
The flux normalizations of the instruments, however, still show 
large deviations among instruments. This is likely caused by 
insufficient modeling of the instrumental deadtimes, the decon- 
volution of the coded mask images, and the calibration of the 
off-axis effective area. 

The high total signal above lOkeV in principle allows us 
to constrain the continuum shape to a higher degree than has 
been possible with RXTE or BeppoSAX alone. The presence 
of an additional, comparatively weak Comptonized component 
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cannot be constrained. In two of the reported cases the ad- 
ditional Comptonized component peaks at ~ 1 ke V, and thus 
only affects the spe ctrum at <10keV JDi Salvo et alJ EoOll 
Fron tera et alJl200lL BeppoSax data). Ginga and OSSE data 
were analyzed using different Comptonization models from 
those applied here an d direct comparisons cannot be made 
dGierlinski et alJ Il997l kTj ~ HOkeV, n ~ 1.7, kT 2 ~ 
45 keV, T2 ~ 6). Int erestingly , the use of a directly comparable 
approach jMaccarone & Connil 120031 eqpair, RXTE/OSSE 
data) also results in a good spectral description using a single 
thermal Comptonization component and parameters fully con- 
sistent with our non-SPI fits. This might indicate a good stabil- 
ity of the hard state spectral shape above lOkeV on a time scale 
of years. 

Thus it is quite clear that apart from model details our 
overall description of the hard state spectrum above 10 keV 
- i.e., with thermal Comptonization and a covering factor 
of the reflector of Q./2n < 0.3 - is principally the same 
as found previously (Dove et a l]ll998t lo ierlinski et al. 1997; 
Pi Salvo et al!200lllFrontera et all200lHMaccarone & Coppil 
2003). Comptonization plus reflection models, featuring a par- 
tial covering of the reflector or a geometrically more compli- 
cated accretion disk configuration, are thus still feasible ex- 
planations of the X-ray spectrum. Direct comparisons between 
these models and other alternative explanations, e.g., involving 
jet physics, still need to be performed. The decrease of the ob- 
served black body flux in the PCA data could be caused by the 
accretion disk fading away after the 2002 soft state. 

We stress that this paper only contains the preliminary 
fits employing the earliest available response models for 
INTEGRAL. These models, however, already show a remark- 
able consistency among independent detectors over a range 
from 1 0-200 keV, with the remaining issues being mainly re- 
lated to the flux calibration of INTEGRAL. Despite these 
limitations we have shown that Comptonization models al- 
ready provide a good description of the joint 3-250 keV 
INTEGRAL-RXTE data. As the calibration of INTEGRAL 
improves, tests using more advanced models than those em- 
ployed here will become possible. Since the INTEGRAL 
data will cover the important range above 250 keV, where 
Klein-Nishina effects become important for the shape of the 
Comptonization continuum and where the non-thermal tails 
start to dominate, it is expected in the near future that these 
data will allow us to constrain the important physical parame- 
ters responsible for the generation of the X-ray spectrum of the 
hard state of Galactic black holes. 
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